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Abstract 

The nonequilibrium thermodynamics feature of a Brownian motor 
operating between two different heat baths is explored as a function 
of time t. Using the Gibbs entropy and Schnakenberg microscopic 
stochastic approach, we find exact closed form expressions for the free 
energy, the rate of entropy production and the rate of entropy flow 
from the system to the outside. We show that when the system is 
out of equilibrium, it constantly produces entropy and at the same 
time extract entropy out of the system. Its entropy production and 
extraction rates decrease in time and saturate to constant value. In 
long time limit, the rate of entropy production balances the rate of 
entropy extraction and at equilibrium both entropy production and 
extraction rates become zero. Furthermore, via the present model, not 
only many thermodynamic theories can be checked but also the wrong 
conclusions that are given due to lack of exact analytic expressions will 
be corrected. 


1 Introduction 

Equilibrium thermodynamics is well established and intensively studied dis¬ 
cipline. However, its application is limited since most systems in nature 
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are far from equilibrium and this warrant the need to advance the nonequi¬ 
librium thermodynamics. On other hand, nonequilibrium thermodynamics 
deals with systems which are nonhomogeneous where the systems thermody¬ 
namic quantities such as entropy, and free energy strictly rely on the system 
parameters in complicated manner. As a result nonequilibrium thermody¬ 
namics is still under investigation. In the last few decades several studies 
have been conducted to explore the nonequilibrium features of systems out 
of equilibrium E :2l El Oh Particularly for system where its dynamics is gov¬ 
erned by a master equation, Boltzmann-Gibbs nonequilibrium entropy along 
with the entropy balance equation serves as a basic tool to understand the 
nonequilibrium thermodynamic features CU21 El- 

Recent study by Schnakenberg also reveals that various thermodynamics 
quantities such as entropy production rate can be rewritten in terms of local 
probability density and transition probability rate [3]. Based on this micro¬ 
scopic stochastic approach, many theoretical studies have been conducted 

m ei ini ei Ei m uni nn usi usi mm These theoretical works confirmed 
that systems which are out of equilibrium constantly produce entropy and at 
the same time extract entropy out of the system. In long time limit, the rate 
of entropy production balances the rate of entropy extraction and at equilib¬ 
rium both entropy production and extraction rates become zero. However, 
most of these previous studies focused on exploring the thermodynamic prop¬ 
erties of the nonequilibrium system at steady state DEI. In this work we 
present a simple model which serves as a basic tool for a better understanding 
of the nonequilibrium statistical physics not only in the regime of nonequi¬ 
librium steady state (NESS) but also at any time t. We explore the short 
time behavior of the system either for isothermal case with load or in general 
for nonisothermal case with or without load. Many mathematical theories 
can be independently checked via the present model [I]. For instance, we 
show that the entropy balance equation always satisfied for any parameter 
choice. Moreover, the first and second laws of thermodynamics are rewritten 
in terms of the model parameters. Several thermodynamic relations are also 
uncovered based on the exact analytic results. 

At this point we want to stress that in this work, we extend (reconsider) 
the previous work nu and uncover far more results. We End exact closed 
form expressions for the free energy dissipation rate F, the rate of entropy 
production e p and the rate of entropy flow from the system to the outside 
h d . The dependence for the change in total entropy, free energy and entropy 
production on model parameters is studied. Since closed form expressions 
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are obtained for all thermodynamic quantities which are under consideration, 
we are able to extract thermodynamic information at any time t. We believe 
that even though, a specific model system is considered, the result obtained 
in this work is generic and advances nonequilibruim thermodynamics. 

Moreover it is found that the entropy S attains a zero value at t = 0; 
it increases with t and then attains an optimum value, ft then decreases as 
t increases further. In the limit t —> oo, S approaches a certain constant. 
Far from equilibrium, in the presence of nonuniform temperature or non-zero 
load, e p and hd decrease with time and approach their steady state value. 
At steady state, e p = hd > 0. At equilibrium, (for isothermal case and zero 
load), e p = hd = 0. Moreover, when the heat exchange via kinetic energy 
is included, we show that even at quasistatic limit e p ^ 0 or E p ^ 0. This 
also implies that since the motor is arranged to undergo a biased random 
work on spatially arranged thermal background, Carnot efficiency or Carnot 
refrigerator is unattainable; there is always an irreversible heat flow via the 
kinetic energy |lT, 18] [T9], 20. 211. 

The rest of paper is organized as follows: in Section II, we present the 
model. In Section III, we derive the expressions the rate of entropy pro¬ 
duction e p and the rate of entropy flow from the system to the outside hd- 
In section IV we explore the free energy F as a function of t. The role of 
particle recrossing on the entropy production rate will be studied in section 
V. Section VI deals with summary and conclusion. 



Figure 1: (Color online) Schematic diagram for a Brownian particle walking 
in a discrete ratchet potential with load. Sites with red circles are coupled 
to the hot reservoir (X),) while sites with blue circles are coupled to the cold 
reservoir (T c ). Site 2 is labeled explicitly and d is the lattice spacing. 
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2 The model 


Consider a single Brownian particle that hops along one dimensional discrete 
ratchet potential with load ||13] 

Ui = E[i(mod) 3 — 1] + ifd (1) 

that coupled with the temperature 

_ , T h , if E[i{mod )3 - 1] = 0; , , 

4 T c , otherwise ; ' ' 


as shown in Fig. 1 . The potential E > 0, / denotes the load and i is an 
integer that runs from — oo to oo. Th and T c denote the temperature for the 
hot and cold reservoirs, respectively. Moreover, the one dimensional lattice 
has spacing d and in one cycle, the particle walks a net displacement of three 
lattice sites. 

The jump probability for the particle to hop from site i to i + 1 is given 
by Ye~ AE ^ ksTi where A E = U l+ \ — Ui and T is the probability attempting 
a jump per unit time, fc# designates the Boltzmann constant and hereafter 
ks, T and d are considered to be a unity. When the particle undergoes a 
biased random walk, it is assumed that first it decides which way to jump 
(backward or forward) with equal probability obeying the metropolis algo¬ 
rism. Accordingly, when A E < 0, the jump the jump definitely takes place 
while A E > 0 the jump takes place with probability exp(—A E/Tf). 

The model exhibits identical behavior with a spin-1 particle system H7]. 
The master equation which governs the system dynamics is given by 


dP n 

dt 


^ ^ ( Pnn'Pn' Pn'nPri) ■ 1 ) 2,3 

n^n' 


(3) 


where P n ’ n is the transition probability rate at which the system, originally 
in state n, makes transition to state n!. Here P n > n is given by the Metropolis 
rule. For instance, P 2 1 = \e~~( E+ P/ Tc , P \2 = |, P 32 = \e~^ E+ P^ Th and 
-P 23 = I- The rate equation for the model can then be expressed as a matrix 
equation ^ = P p where p = (pi,P2iP3) T ■ P is a 3 by 3 matrix which is given 
by 



( 4 ) 
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as long as 0 < / < 2 E/d. Here p = e -E / Tc , v = e~ E ! Th , a = e _ ^ d//Tc 
and b = e~^ d ^ Th . Note that the sum of each column of the matrix P is 
zero, P nm = 0 which reveals that the total probability is conserved: 
(d/dt) J2 n Pn = d/dt( 1 T ■ p) = 1 T • (P p) = 0. 

For the particle which is initially situated at site i = 2, by solving Eq. 
(3), we find the time dependent normalized probability distributions for pi(t), 
p 2 (t) and p 3 (t) as shown in Appendix A. The expression for the velocity V{t) 
at any time t as well as the rate of heat flow from the hot reservoir Qh(t) and 
the rate of heat flow into the cold reservoir Q c (t ) are also given in Appendix 
A. 

Hereafter all figures are plotted by taking dimensionless quantities e = 
E/T c , A = fd/T c and r = ^ — 1. We also introduce dimensionless time 
t = Tt and after this the bar will be dropped. 


3 Derivation of Entropy and entropy produc¬ 
tion 

The relation between the internal energy U(t), entropy S(t) and free energy 
F(t) is well known for isothermal system (uniform temperature T) which 
satisfies detailed balance condition, and can be written as 

U\Pi(t)} = 


= 


and the corresponding fundamental entropy balance equation is given by 

S{t) = e p (t)-h d (t) (7) 

where S is the Gibbs entropy given by 

N 

S\pi(t)] = -^Pilnpi. (8) 


-Th d {t) = -H d (t ) 


(5) 


-T XX diPji - PjPij) ^ ^ , 

i>j k' './ / 


/V 


- Te p (t ) = —E P (t) 

-T piPji - PjPij) In 


( 6 ) 


i>j 


P> 1> , 


3 l 


PjPi 
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Here hdit) and e p (t) designate the term which is related to the heat dissi¬ 
pation rate and the instantaneous entropy production (e p > 0), respectively. 
The entropy balance equation (7) can be rewritten as S T (t ) = E p (t) — Hdit ) 
where S T (t) = TdS{t)/dt. Recently for isothermal system which is driven 
out of equilibrium, the above relations have been extended via phonological 
approach [T] . 

Next, we examine whether the well-known thermodynamic relations, which 
are valid for equilibrium system, are still obeyed for the model system which 
is driven out of equilibrium due to inhomogeneous thermal arrangement 
Th 7 ^ T c or non-zero external load / ^ 0. For the three state model we 
present, 

3 

S\pi(t)] =-^pilnpi, (9) 

2=1 


and 


h d 


^ ] iPi.Pji 

i>j 


PjPij) ln 



(-P1-P21 +^2^12) In 
(-P2-P32 + P 3 -P 23 ) In 


PvA 

P21) 

P 2 A 

P 32 J 


(P3-P13 - P1P31) in 



+ 

+ 


^(PiPn ~ Pj p v ) ln ( 77 ^) 

i> j \Pj r HJ 

(-P1P2I + P2P12) In + 

(~pA + P 3 P2 3 )h,(^PPj + 

(P3P13 - P1P31) in (~ 7 =n) 


( 10 ) 


(ii) 


where the indexes i = 1...3 and j = 1...3. ft is important to note that via the 
expressions pi(t), P 2 (t) and p 3 (f) that are shown in Appendix A and using 
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the rates (see Eq. (4)), 

P21 = \e - (B+/,/T - Pi2 = i P:,2 = fi^E+n/n 

As = i Pis = 5, P 31 = i e - (2E - fl/T - ( 12 ) 

the thermodynamic quantities which are under investigation can be evalu¬ 
ated. 

For our nonequilibrium system where T h T c and f ^ 0, regardless of 
our parameter choice, we find that the fundamental entropy balance equation 


S(t) = e p (t)-h d (t) 


(13) 


is always satisfied at any time t. Moreover h d {t) is found to satisfy the 
relation 


h d {t) 


-Qhjt) Q c (t) 

T h T c 


(14) 


Next we study how S, S, and e p and h d vary in time. 

Entropy. — The entropy of the system exhibits an intriguing parameter 
dependence. Exploiting Eq. (9) one can see that (for any parameter choice) 
when t —> 0, S —> 0. As time increases S increases and attains an optimum 
value. As t further increases, S decreases. In the limit t — * 00 , the total 
entropy approach its steady state value S — > — Y^=iPi l n pf- Particularly, 
for the case where / = 0, one can get a simplified expression by expanding 
S in small t regime and after some algebra one gets 

1 -Up Up Up 

S = -e T h (1 + e T h + e T h ln[2])t (15) 

-up 

-Up £ T h £0 

~2 e Th - ln M - eTh ln MK 

Eq. (15) clearly depicts that in the limit t —> 0, S —> 0. In Fig. 2a, we plot 
the entropy S as a function of t for parameter choice of Uq = 2.0 , r = 2 and 
A = 0 (red line) and A = 0.5 (black line). The figure once again shows that 
the entropy attains an optimum value at particular time t. 

One can note that the system sustains a non-vanishing current (the sys¬ 
tem is out of equilibrium) as long as T\ T c and Uq 0 both in the presence 
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or in the absence of load. For isothermal case, the system is driven out of 
equilibrium only in the presence of load / > 0. This implies that the sys¬ 
tem relaxes to equilibrium in the absence of load and when Th = T c . In 
the absence of both external load and bistable potential Uq = 0, the particle 
undergoes a random walk on lattice. In long time limit, the system may 
approach equilibrium as long as the heat exchange via kinetic energy is ne¬ 
glected (even if a distinct temperature difference is retained between the hot 
and cold reservoirs). Thus when / = 0 and in the limit Uq —!> 0 (approaching 
equilibrium), Eq. (9) converges to 

S{t) = \e^(- 2(— l + e f)ln(l-^)) 

1 -31 , . 3t . . .1 2 -31 . . , 

--e 2 ((2 + e 2 ) hi (- + -e 2 )). (16) 

As it can be seen from Eq. (16), in the limit S —> 0 as t — y 0 and when 
t —> oo, S —> In[3]. Note that at stationary state, S is constant implies that 
AS = 0. 



Figure 2: (Color online) (a) Total entropy S(t) versus t for a given e = 2, 
r = 2.0 and A = 0.(b) Plot of S(t) (dotted green line), e p (t) (dashed red line) 
and hd(t) (black solid line) as a function of t for a given e = 2, r = 2.0 and 
A = 0. 

Entropy production rate .— Let us focus on the rate of entropy production 
e p , the rate of total entropy S(t) and the rate of entropy flow from the system 
to the outside hd(t). The plot of S(t) (dotted green line), e p (t) (dashed red 
line) and hd(t) (solid black line) as a function of t is depicted in Fig. 2b 
for a given values of e = 2, r = 2.0 and A = 0. The figure indicates that 
e p and h d (t) decrease to their steady state values e p = h d > 0.04 as time 









0.010 



A 

Figure 3: (Color online) (a) Plot of S(t) (green line), e p (t) (dashed red line) 
and h d (t) (solid line) as a function of t for a given e = 2, r = 2.0 and t = 10 6 
(steady state). Here e p (t) = h d (t). 

progress. Exploiting Eqs. (10), (11) and (13), one can also see that far from 
equilibrium e p > 0 and h d > 0. When t —> 0, e p becomes much greater than 
h d {t) and as time increases, in a certain time interval, h d {t) > e p . In general, 
as time progresses e p and h d decrease and approach their steady state value. 
At steady state, e p = h d > 0. At equilibrium, for isothermal case and zero 
load, e p = h d = 0. Particularly expanding e p and h d {t) in the small time 
regime, we find 

1 zEsl 2 zEsl e ~£a 

= 2 e Th - e Th “ ( eTh ) ln W) ( 17 ) 

and 

— 1 -Vfi 1 ~ u 0 Vq 

hd = — ln[e Tc ] - -e T h ln[e T fe]. (18) 

On the other hand in the limit Uq —> 0 and / = 0 (approaching equilibrium), 
hd = 0 for any t while 

. — 3 1 . — 1 -(■ 6 2 , . . 

e P = —e 2 ln[—-(19) 

2 + e 2 

Here e p > 0 for small t and decreases (the system relaxes to equilibrium) as 
time increases. In the limit t —> oo, e p = h d = 0. 

The quasistatic limit of the system corresponds to the case where the 
velocity approaches zero V(t) = 0 in the vicinity of the stall force or in the 
absence of load when Uq —» 0. The general expression for the stall force 
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is lengthy and will not be presented. However, at steady state, the stall 
force reduces to / = E(j± — 1)/(2^ + 1). At quasistatic limit, regardless of 
any parameter choice, we find e p = hdit) = 0. This can be appreciated by 
plotting e p or hdit) as a function of load in long time limit. Figure 3 exhibits 
that both e p and hdit) decrease as the load increases and attain a zero value 
at stall force A = 0.4. As the load further increases, e p and hdit) step up. On 
the contrary, 5 decreases with time. At stall force, 5 = 0. 



Figure 4: (Color online) (a) Plot of Ae p = e p (t) — e p (0) (red line) and 
A hd = hdit) — hdit)) (Black line) as a function of t for a given e = 2, r = 2.0 
and A = 0. (b) A5(f) — S(t) — 5(0) versus t for a given e = 2, r = 2.0 and 
A = 0. 


Since closed form expressions for 5(f), e p (t) and hdit) are obtained for any 
time f, the analytic expressions for the change in entropy production, heat 
dissipation and the change in the total entropy can be found analytically via 


and 



( 20 ) 


( 21 ) 



3 3 

= -^2pi(t)lnPi(t)+ ^2pi(0)lnpi(0) (22) 

2=1 2=1 
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where AS = Ae p — Ah d and the indexes i = 1...3 and j = 1...3. The 
expressions for A h d , AS and Ae p are lengthy. 

Here care must be taken since e p (t) > 0 or hdit) > 0 does not imply that 
Ah d = h d (t)-h d (t 0 ) > 0, AS = S(t)-S(t 0 ) > 0 or Ae p = e p (t)-e p (t 0 ) > 0. 
Rather since hd, S and e p are decreasing function of t, Ah d , AS and Ae p can 
take even negative values depending on the interval between £ 0 and t. If the 
system reaches at steady state or stationary state and if the change in these 
parameters are taken at the steady state regime, then A h d = 0, AS = 0 or 
Ae p = 0. However this does not indicate that second law of thermodynamics 
is violated. 

In fact the second law of thermodynamics is always preserved as it can 
be seen by exploiting the exact analytic expressions (Eqs. (20-22). In reality, 
once the motor starts operating, entropy will be accumulated in the system 
starting from t — 0 and as time progresses, more entropy will be stored in 
the system even though some entropy is extracted out of the system. Hence 
if the change in these parameters is taken between t — 0 and any time t, 
always the inequality A h d = h d (t) — h d ( 0) > 0, AS = S(t ) — 5(0) > 0 or 
Ae p = e p (t) — e p (0) > 0 holds true and as time progresses the change in this 
parameters increases as shown in Fig. 4a. As a matter of fact, in small t 
regimes, e p becomes much greater than hdit) (see Fig. 2b) revealing that the 
entropy production is higher (than entropy extraction) in the first few period 
of times. As time increases, more entropy will be extracted h d (t) > e p . Over 
all, since the system produces enormous amount of entropy in initial times, 
in latter time or any time t, Ae p > Ah d and hence AS > 0 (see Fig. 4). 


4 Free energy dissipation rate and Free en¬ 
ergy 

In order to relate the free energy dissipation rate with E p {t) and Hdit) let us 
now introduce H d {t) for the model system we considered. The heat dissipa¬ 
tion rate is given by H d (t) = — Qh(t ) + Q c [t). Using Eqs. (56) and (57), one 
finds 


H d {t) 


— Qh{t ) + Qc(t) 

T o(Pi P H - Pj p ii) ln 

i>j 



(23) 
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Equation (23) is notably different from Eq. (10) due to the term Tj. Similar 
relation has been used by Hao. et. al for the isothermal case |T]. We rewrite 
Eq. (23) as 


where 


and 


H d {t) = ^2 T i _ Pj p ij) ln 


i>j 


-pjPij) In 


i>j 


P 

1 31 

P 

1 13 

Pi P 3 


31 


J2 _ Pj P ij) 111 


i>j 


PjPij 


Pi 


Pi 


= Ep(t) — S T (t) 


Ep(t) = ^TjiPiPji-PjPij) In 


i>j 


PiPj 


3i 


Pj P i 
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S T (t) = V 7}(/>,/ , ; , - PjPij) In 


i>3 


(24) 

(25) 

(26) 


Here the indexes i = 1...3 and j = 1...3. Now we have entropy balance equa¬ 
tion S T (t) = E p (t ) — Hd(t) for our model system. Note that for isothermal 
case Eqs. (24), (25) and (26) converge to Thd(t ), Te p (t ) and S T (t) = TS(t ), 
respectively. All the above analysis indicates that E p (t), Hd(t ) and S T (t ) con¬ 
tain a term ( P 2 P 32 — P 3 P 23 ) which is associated with the rate of heat taken 
out of the hot reservoir and two terms (P 2 P 12 ~ P 1 P 21 ) an d (P 3 P 13 — P 1 P 31 ) 
which are associated with the rate of heat given to the cold reservoir. 

One notable difference between e p (t) / hd(t) and E p (t) / Hdif) is that e p (t) / hd(t) 
always approaches unity in long time limit while E p (t)/Hd(t ) approaches 
unity only for isothermal cases with or without load as shown in Fig. 5. 
Figure 5a depicts the plot of e p (t)/hd(t) as a function of t for a given e = 2, 

/ = 0.8 r = 2.0 (red line) and r = 1.0 (black line). Figure 5b shows the Plot 
of E p (t)/Hd(t) as a function of t for a given e = 2, / = 0.8 and r = 2.0 (red 
line) and r = 1.0 (black line). 
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Figure 5: (Color online) (a) Plot of e p (t)/h d (t ) as a function of t for a given 
e = 2, / = 0.8 r = 2.0 (red line) and r = 1.0 (black line), (b) (a) Plot of 
E p (t)/H d (t) as a function of t for a given e = 2, / = 0.8 r = 2.0 (red line) 
and r = 1.0 (black line). 


The second law of thermodynamics can be written A S T = A E p — A H d 
where A S T , A E p and A H d are very lengthy expressions which can be eval¬ 
uated via 


and 


AH, = 


A E p = 


AS t = 


[ ['52 T j(Pi p ji-Pj p ij) in 

■ h ° \i>j 

[ ( ln 

Jt ° \i>j 



PiPji 
PjPij - 


dt 


f ('E T ite p ii-p i Pii) 

Jt 0 V; ( VPj, 


dt. 


(27) 


(28) 


(29) 


In Fig. 6a, we plot A E p = E p (t) — E p { 0) (red line), A H d {t) = H d (t) — ^(0) 
(green line) and A S T (t) = S T (t ) — S' T (0) (black line) as a function of t for 
given e = 2, / = 0.8 and r = 2.0. The figure exhibits that as long as t > 0, 
A E p > A H d and hence A S T > 0. 

On the other hand, the total internal energy U(t) is the sum of the internal 
energies 


U\pi(t)] = 


2=1 


= Pl (t)(-E)+p 3 (t)(E) 


(30) 
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Figure 6: (Color online) (a) Plot of A E p (t (red line) A H c i(t) (green line) 
and A S T (t) (black line) as a function of t for a given e = 2, / = 0.8 and 
r = 2.0 (red line), (b) Plot of U(t) as a function of t for a given the same 
parameter choice. 


while the change in the internal energy is given by 

AU = U\pi(t)} — U\pi(0)} 

= E(Ps(t) ~p 3 (0) +Pi(0) -pi(t)). (31) 

As the particle walks along the reaction coordinate, it receives some heat 
from the hot reservoir and gives part of it to the cold bath. The remaining 
heat will be spent to do some work W(t) = fV(t ) against the external load. 
Hence we verify the first law of thermodynamics 

U[Pi(t)] = -^(PiPji-PjPi^iui-Uj) 

i>j 

= Qh{t)-Q c {t)-W{t ) 

= ~(H d (t) + fV(t)). (32) 

The change in the internal energy (Eq. (31)) can be also rewritten as 

AU = - J (Q h (t)-Q c (t)-W{t)}dt 

= J' (-(H d (t) + fV(t)j)dt 

= A Q h - A Q c - AW. (33) 

Note that similar relation is derived in the work [23] for the isothermal case. 
Fig. 6b depicts the plot of U{t) as a function of t for the same parameter 
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Figure 7: (Color online) (a) The free energy dissipation rate F versus t for 
a given e = 2, r = 2.0 and A = 0.8. (b) Plot of S T (t) (dotted line), E p (t) 
(dashed line) and H^it) (solid line) as a function of / for a given e = 2, 
r = 2.0 and A = 0. 


choice of e = 2, / = 0.8 and r = 2.0. The figure shows that as t —> oo, 
U(t) —)■ 0 as expected. 

Next let us find the expression for the free energy dissipation rate F. For 
the isothermal case, the free energy is given by F = U — TS and we can still 
adapt this relationship to nonisothermal case to write 

F(t) = U-S T (t). (34) 

Substituting Eqs. (23) and (25) in Eq. (32) leads to 

m + Ep(t) = U(t) + H d (t) = —fV(t) (35) 

which is the second law of thermodynamics. Note that in the absence of load, 
U(t) = —H d (t) and consequently E p (t) = —F(t ). The dependence of the free 
energy on time t can be explored by exploiting Eq. (35). In general F < 0 
and approaches to zero in the long time limit (see Fig. 7a). In order to get 
a very simplified expression, let us expand Eq. (35) in small time regime for 
the case where / = 0 to get 

-Uo 

• 1 — ^0 t/p f T c 

F = -e"ET (~T h In [2] + eE ln[-y ]T C ) + 

1 zEo £o 

-e T h (e T h T c ln[t] + T h ln[t]). (36) 

When t —> 0, F diverges. On the other hand when U$ —> 0 and / = 0 


15 








(approaching equilibrium), we get 

j’ = e^ i (T c + T»)ln[i^i]. (37) 

2 + e^ 

As it can been seen in Eq. (37), in the limit t —y oo, F —> 0 and when t —> 0, 
F diverges. 

The change in the free energy also is given by 

A.F = -jf (fV(t) + E p {tj)dt (38) 

= f (u(t) + H d (t) - 4(t)) dt 
= AU + AH d - AE P . 

At quasistatic limit where the velocity approaches zero V{t) = 0, E p (t) = 
0 and H d {t) = 0 (see Fig.(7b)) and far from quasistatic limit E p > 0 which 
is expected as the engine operates irreversibly. Far from stall force, E p (t ) 7 ^ 
H d (t) as long as a distinct temperature difference between the hot reservoirs 
is retained. However, for isothermal case, at steady state, E p (t) = Hd(t) 
which implies A F = AU. 

Finally, in the work [lj, the second law of thermodynamics is rewritten 
in terms of the housekeeping heat and excess heat by considering isothermal 
case. For the model system we consider, when the particle undergoes a 
cyclic motion, at least it has to get fV (t) amount of energy rate from the 
hot reservoir in order to keep the system at steady state. Hence fV(t) is 
equivalent to the housekeeping heat Qhk and we can rewrite Eq. (35) as 

F(t) + E p (t) = U(t) + H d (t) = -fV(t ) = -Q hk (39) 

while the excess heat Q ex part is given by 

Qex = H d — Q hk . (40) 

For isothermal case, we can rewrite the second law of thermodynamics as 
S T = E P -H d =-F-Q ex (41) 

and 

F = Qhk — E p . (42) 


16 



5 Irreversibility due to particle recrossing along 
the boundaries 

In the previous section we have seen that at quasistatic limit, the system is 
reversible and hence e p = 0 or E p = 0 even if T h y T c . However so far the 
rate of heat loss due to particle recrossing at the boundary between the hot 
and cold reservoirs is not included. If the heat exchange via kinetic energy 
is included, even at quasistatic limit e p ^ 0 or E p ^ 0. This also implies 
that since the motor is arranged to undergo a biased random work on spa¬ 
tially arranged thermal background, Carnot efficiency or Carnot refrigerator 
is unattainable; there is always an irreversible heat flow via the kinetic energy 

pn na m, he m\. 



0.0 0.2 0.4 0.6 0.8 1.0 
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Figure 8: (Color online) Plot of E p {t) (dashed black line) and E*(t ) (red 
line) as a function of A for a given e = 2, r = 2.0 and t = 10 6 . 

When the particle by chance jumps from the cold to hot reservoir, it 
absorbs, ks{Th — T c )/2 amount of heat from the hot bath. When it hops 
back to the cold reservoir, it dumps back this ks{Th — T c )/2 amount of heat 
to the cold bath. This irreversible heat flow is always one way, i.e.; from the 
hot to the cold baths and it cannot be recovered. Hence we can conclude 
this heat transfer from the hot to cold reservoirs depends on how often the 
particle jumps from the cold to hot heat baths and it can be written as 

Qirr = (^1^21 + £>3-^23) {Th ~ T c )/ 2. (43) 

The heat exchange via kinetic energy does not affect since the whole heat 
taken from the hot reservoirs goes to the cold reservoir. This also implies 
that the whole heat dumped to the cold reservoirs contributes to the internal 
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entropy production and hence for any parameter choice E p > 0 as long as 
T h ^ T c . The heat loss due to particle recrossing only contributes to the 
internal entropy production and we infer the new entropy production rate to 
be 


E* — Ep + Qirr . (44) 

One can then rewrite the thermodynamic relations that derived in the previ¬ 
ous section in terms of E*. Here E* = 0 only when Th = T c . E* as a function 
of the load A is plotted in Fig. 8. The figure depicts that even at quasistatic 
limit e* ^ 0. The same figure shows that if the heat exchange via kinetic 
energy is omitted, E p = 0 at stall force A = 0.4. Note that when Uo and 
/ = 0, the particle undergoes unbiased random walk problem. Its average 
velocity is zero but its speed is nonzero. If the particle by chance hops from 
the cold to hot reservoirs, it absorbs an energy and latter dumps this energy 
to cold bath. So even in this case the system is far from equilibrium. 

6 summary and conclusion 

In this work we present a paradigmatic model which serves a basic tool for a 
better understanding of the nonequilibrium statistical physics not only in the 
regime of nonequilibrium steady state (NESS) but also at any time t. Not 
only the long time behavior of the system is explored but also we investigate 
the short time behavior of the system either for isothermal case with load or 
in general for nonisothermal case with or without load. 

Based on the Gibbs entropy and Schnakenberg microscopic stochastic ap¬ 
proach, we get exact analytic expressions for the thermodynamic quantities 
that are under investigations. We show that whenever the system is out 
of equilibrium, it constantly produce entropy and at the same time extract 
entropy out of the system. At steady state, the rate of entropy produc¬ 
tion balances the rate of entropy extraction and at equilibrium both entropy 
production and extraction rates become zero. Moreover we show that the 
entropy balance equation [22] always satisfied for any parameter choice. Fur¬ 
thermore, the first and second laws of thermodynamics are rewritten in terms 
of the model parameters. Several thermodynamic relations are also uncov¬ 
ered based on the exact analytic results. 

The exact analytic expressions are also obtained for the free energy dissi¬ 
pation rate F, the rate of entropy production e p and the rate of entropy flow 
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from the system to the outside h d . Particularly, the analytical study reveals 
that the entropy S attains a zero value at t — 0; it increases with t and then 
attains an optimum value. It then decreases as t increases further. In the 
limit t —>• oo, S approaches a certain constant. Far from equilibrium, in the 
presence of nonuniform temperature or non-zero load, e p and h d decrease with 
time and approach their steady state value. At steady state, e p = h d > 0. 
At equilibrium, for isothermal case and zero load, e p = h d = 0. Moreover, 
when the heat exchange via kinetic energy is included, we show that even at 
quasistatic limit e p ^ 0 or E p ^ 0. This also implies that since the motor 
is arranged to undergo a biased random work on spatially arranged thermal 
background, Carnot efficiency or Carnot refrigerator is unattainable; there is 
always an irreversible heat flow via the kinetic energy. 

In conclusion, in this work we present a simple model which serves as a 
basic tool for better understanding of the nonequilibrium statistical physics 
not only in the regime of nonequilibrium steady state (NESS) but also at 
any time t. The present model also serves as a tool to check many elegant 
thermodynamic theories. Based on this exactly solvable model, we uncovered 
several thermodynamic relations. We believe that even though, a specific 
model system is considered, the result obtained in this work is generic and 
advances nonequilibruim thermodynamics. 
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Appendix A 

In this Appendix we will give the expressions for pi(t), p 2 (t) and p${t) as well 
as V(t), Qh(t) and Q c (t) For the particle which is initially situated at site 
i = 2, the time dependent normalized probability distributions after solving 
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Eq. (3) are 


where 


Pi(t) = Ci 


a( 2 + vb) 


p (p + (a 2 + p) vb) 


+ 


(a+a^-j-^ ) t 


c 2 e 


-1 + 


a(—1 + a/x)\ 

—/i 2 + az/6 J ’ 


P2(t) = 


_ lt(-2-vb) „ a e 2 “ ( _1 + a P) 


Cl 


-c 3 e 2 v ' - c 2 - 
(2a 2 + p) 


—p 2 + az/6 


Pzit) = ci + c 2 e 


p + (a 2 + /i) z/6’ 

(" a+a^ p+p^^t -i 

1 + c 3 e~^- 2 - vb) 


Cl = 

C2 = 

C 3 = 


+ 


p(p + (a 2 + p) vb) 

(a + a 2 p + p 2 ) (2 + vb) ’ 
a 

(a + a 2 p + p 2 ) (-1 + ) 

p (p + a 2 z/6 + pvb) 

(a + a 2 p + yU 2 ) (2 + vb) ~*~ 

a 

(a + a 2 fi + M 2 ) (-1 + “IJ+aHt ) 


(45) 

(46) 

(47) 

(48) 

(49) 

(50) 


The J2i=i Pi(t) = 1 revealing the probability distribution is normalized. In 
the limit of t —» oo, we recapture the steady state probability distributions 


Pi 

Pi 

Pi 

shown in the work [15] . 


a 

a + a 2 p + p 2 ’ 

p (2 a 2 + p) 

(<a + a 2 p + p 2 ) (2 + bv) ’ 
p (p + b (a 2 + /i) y) 

(a + a 2 /i + p 2 ) (2 + bv) 


(51) 
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(53) 


20 


















The velocity V(t) at any time t is the difference between the forward 
V+it) and backward V~{t ) velocities at each site i 

V(t) = £(v; + (t) - v;-(i)) (54) 


1 = 1 


= (piAl - P2P12) + (P2P32 ~ P3P23) 
(P.'i P\ 3 -Pl^il)- 


Exploiting Eq. (14), one can see that the particle retains a unidirectional 
current when / = 0 and Th > T c . For isothermal case Th = T c , the sys¬ 
tem sustains a non-zero velocity in the presence of load / 7 ^ 0 as expected. 
Moreover, when t —>■ 00 , the velocity V(t) increases with t and approaches 
to steady state velocity 




b ( bau ~ g) 

2(2 + ub) ^1 + a/i + 


(55) 


The heat per unit time taken from the hot reservoir has a form 
Qh(t) = (E + /)(p 2 P 32 ~ P 3 P 23 ) 

= T h (p 2 P 32 -P3^23)ln( : ^). (56) 

2 23 

On the other hand, the heat per unit time given to cold reservoir is given by 

Qcif) = (E + /)(p 2 Pl 2 — P1P21) + 

(2 E - f)(p 3 Pi 3 -P 1 P 31 ) 

= T c (p 2 P 12 - Pl P 21 ) ln(^) + 

P 21 

Tc(P3Pi3-PiP3i)H^) (57) 

Oil 

As t —> 00 , both Qh(t) and Qc(t) evolve in time to their corresponding steady 
state values 


Qh = (E + f) 


/i (bav - P) 


2(2 + vb) ^1 + ap + 


and 


Qc = [E-VY 


T {bav - E) 


2(2 + vb) (l + ap+E 


(58) 


(59) 
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